PENTAQUARK SEARCH AND OTHER MULTIQUARK CANDIDATES AT 

BES 



XIAOYAN SHEN (FOR THE BES COLLABORATION) 
Institute of High Energy Physics, Chinese Academy of Sciences, P. R. China 

Results are presented on ip{2S) and J/ip hadronic decays to KgpK~n and K^pK+n final states from 
data samples of 14 million ip(2S) and 58 million J/ip events accumulated at the BES II detector. No 
0(1540) signal, the pentaquark candidate, is observed, and the upper limits are set. We also present 
a narrow enhancement near 2m p in the invariant mass spectrum of pp pairs from J/tp — * -ypp decays, 



an enhancement near the m p 
threshold from J/ip — » pK~A ■ 



- Ma mass threshold and an enhancement near the n%K 

- c.c. decays, based on 58 million J/ip events. 



m,A mass 



1 Pentaquark search 

In recent years, there have been many ex- 
periments reporting the observation of a new 
narrow S = +1 state, the 0(1540), and 
also many experiments claiming the null ob- 
servation of the 0(1540). Therefore, fur- 
ther experimental confirmation is required 
before claiming solid evidence for a S = +1 
baryon resonance. Compared with the above 
experiments, the data accumulated at the 
e + e~ collision experiment BES are relatively 
clean and have less background; therefore it 
is meaningful to investigate the pentaquark 
state O with the hadronic decays of ip(2S) 
and J/ip 4 . 

The pentaquark state 0(1540) in tp(2S) 
and J/ip decays to KgpK~n and KgpK + n 
final states with Kg decaying to tt + -k~ is 
searched for using 14 million ip(2S) and 
58 million J/ip events taken with the up- 
graded Beijing Spectrometer (BES II) 1 lo- 
cated at the Beijing Electron Positron Col- 
lider (BEPC). 

Fig. 1 (left) and Fig. 1 (right) show 
the scatter plots of K~n (Kgp) versus 
K s p (K+n) for ip(2S) and J/ip -> K° s pK~n 
+ KgpK + n modes. No clear O signal is ob- 
served in both ip(2S) and J/ip decays. The 
upper limits are set at the 90% confidence 
level (C.L.) as: 

B(iJj(2S) -> OO -» K^pK-fi+c.c. < 0.84x10" 



□ 



M<K s p(K (GeV/c") 
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B(ip(2S) -> QK-n -> K° s pK-n) < 1.0x10 



-5 



Figure 1. Left: Scatter plot of K n (K®p) versus 
K° s p (K+n) for tp{2S) -» K s pK~n + K°pK+n 
modes. Right: Scatter plot of K~n (K^p) ver- 
sus K%p {K+n) for J/ip -» K%pK~n + K%pK+n 
modes. 



B(ip(2S) -» QK+n -> K s pK + n) < 2.6xl0~ 5 
B{ip(2S) -> K° s pO -> K s pK-n) < 0.60xl0~ 5 
B(<ip(2S) -> K s pO -> KgpK+n) < 0.70 xl0~ 5 . 
B(J/ip -> 00 -> K sP K~n+c.c.) < 1.1 xl0~ 5 
B(J/ip -► QK~n -► K%pK~n) < 2.1 x 10~ 5 
B(J/ip -> @K + n -> K s pK + n) < 5.6 x 10~ 5 
B(J/ip -> K s pS -» KgpK-fi) < 1.1 x 10~ 5 
B( J/ip -> K s p@ -> K s pK + ii) < 1.6 x 10~ 5 . 

2 Near pp threshold enhancement 

in J/ip — > jpp 

5 There is an accumulation of evidence for 
anomalous behavior in the proton-antiproton 
(pp) system very near the M p p — 2m p mass 
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Figure 2. The pp invariant mass distribution 

J/tfi — > ■ypp decays. 



threshold 2 . Based on a sample of 58 million 
J/ip events accumulated at BES, we report a 
study of the low mass pp pairs produced via 
radiative decays 5 . 

The events with a high energy gamma 
ray and two oppositely charged tracks, each 
of which is well fitted to a helix originat- 
ing near the interaction point and within 
the polar angle region |cos#| < 0.8, are se- 
lected. The dE/dx information is used for 
the particle identification and both charged 
tracks should be identified as proton and anti- 
proton. The candidate events are subjected 
to four-constraint kinematic fits to the hy- 
potheses J/ip — > jpp and CLjpp > 0.05 is 
required. 

Fig. 2 shows the pp invariant mass dis- 
tribution for surviving events. Except for a 
peak near M pp — 2.98 GeV/c 2 that is consis- 
tent in mass, width, and yield with expecta- 
tions for J/ip — > 7?7 C , rj c — > pp 3 and a broad 
enhancement around M pp ~ 2.2 GeV/c 2 , 



there is a narrow, low-mass peak near the pp 
mass threshold. 

The low mass region of the pp distribu- 
tion is fitted by an acceptance-weighted S- 
wave Breit-Wigner (BW) function to repre- 
sent the low- mass enhancement and /bk g (<5) 
to represent the background. The mass and 
width of the BW signal function are allowed 
to vary and the shape parameters of fbkg(8) 
are fixed at the values derived from the fit 
to the TT°pp phase-space MC sample. This fit 
yields 928±57 events in the BW function with 



a peak mass of M = 1859 
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MeV/c 



and a full width of V < 30 McV/c 2 at a 
90% confidence level (CL). Here the first er- 
ror in the mass is the statistical error and 
the second is the systematic error which in- 
cludes changes observed in the fitted values 
for fits with different bin sizes, with back- 
ground shape parameters left as free parame- 
ters, different shapes for the acceptance vari- 
ation, different resolutions, and the range 
of differences between input and output val- 
ues seen in the MC study. Using a Monte- 
Carlo determined acceptance of 23%, we 
determine a product of branching fractions 
B{J/ip -► 7X(1859))£pf(1859) -► pp) = 
(7.0±0.4(stat)+Q g(syst)) x 10~ 5 , where the 
systematic error includes uncertainties in the 
acceptance (10%), the total number of J/ip 
decays in the data sample (5%), and the ef- 
fects of changing the various inputs to the fit 

f+24%\ 

I -2%)- 



3 The anomalous enhancements 
near the m p + M\ and vhk + M\ 
mass thresholds in 

J/ip —> pK~A + c.c. decays 

The results from J/ip — > pK~A + c.c de- 
cays are presented 6 . The J/ip —> pK~A 
candidate events are required to have four 
charged tracks with net charge, each of which 
is well fitted to a helix within the polar an- 
gle region | cos#| < 0.8 and with a transverse 
momentum larger than 50 MeV. The com- 
bined TOF and dE/dx information is used 
for the particle identification. Events where 
the p, K~, p and tt + tracks are all unam- 
biguously identified are subjected to a four- 
constraint (4C) kinematic fit with the cor- 
responding mass assignments for each track. 
For events with ambiguous particle identifica- 
tion, all possible 4C combinations are formed, 
and the combination with the smallest x 2 
is chosen. The final x 2 is required to be 
less than 20. The background events from 
J/ip — > pK~Y, are suppressed by requiring 



For Publisher's use 



6 = E miss +l.39M pK miss < 1.69 GeV, where 
Emiss denotes the difference between the 
center-of-mass energy (3.097 GeV) and the 
total energy of the four charged tracks, and 
M p k miss denotes the mass recoiling against 
the proton-kaon system. The monte-Carlo 
studies indicate that the background in the 
selected event sample is at the 1 ~ 2% level 
after above criteria. 

The pA invariant mass spectrum for the 
selected events is shown in Fig. 3(a), where 
an enhancement is evident near the mass 
threshold. The pK~A Dalitz plot is shown in 
Fig. 3(b). In addition to bands for the well 
established A*(1520) and A* (1690), there is 
a significant N* band near the K~A mass 
threshold, and a pA mass enhancement, iso- 
lated from the A* and TV* bands, in the right- 
upper part of the Dalitz plot. 

This enhancement can be fitted with 
an acceptance weighted S-wave Breit-Wigner 
together with a function fps{8) describing 
the phase space contribution, as shown in 
Fig. 3(c), where f PS (8) = N(8 X / 2 + M 3 / 2 + 
«2^ 5 ^ 2 ), 8 — m pA — Trip — m A , and the param- 
eters ai and a,2 are determined from a fit to 
the pK~A MC sample events generated with 
a uniform phase-space distribution. The fit 
gives a peak mass of m = 2075 ±12 MeV, 
a width r — 90 ± 35 MeV and a branching 
ratio 



BR{J/%I> -► R-X)BR(X 
N res /(2eBR(A -> pn)) 



N 



P A) 

= (5.9 ± 1.4) x 10" 



Performing the same analysis on the ip' 
data sample, an evidence of a similar en- 
hancement is observed in ip — > pK~A, shown 
in Fig. 4(a) and (b). If the V(2075) param- 
eters are fixed to the values obtained from 
the J/ip data, i.e., M x = 2075 MeV and 
Tx = 90 MeV, the fit to the ip' data sample 
shows that the enhancement in ip' data devi- 
ates from the shape of the phase space contri- 
bution with a statistical significance of about 
4.0cr, where the significance is estimated from 
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Figure 3. (a) The points with error bars indicate the 
measured pA mass spectrum; the shaded histogram 
indicates phase space MC events (arbitrary normal- 
ization), (b) The Dalitz plot for the selected event 
sample, (c) A fit (solid line) to the data. The dot- 
ted curve indicates the Breit-Wigner signal and the 
dashed curve the phase space Background', (d) The 
cos Bp distribution under the enhancement, the points 
are data and the histogram is the MC (normalized to 
data) 



a comparison of log-likelihood values of the 
fits with and without the X(2075) signal func- 
tion. 

As mentioned above, the pK~A Dalitz 
plot, Fig. 3(b), shows a significant N* band 
near the K~A mass threshold. This band 
corresponds to an enhancement near the 
K~ A mass threshold in the one dimension 
projection of m K - A , shown in Fig. 5(a). The 
5 m K-A~ m K- — toa distribution after the effi- 




Figure 4. Results for t/} 1 



(Gen 

pK~K events: (a) A 



fit (solid line) to the data sample (histogram); the 
dashed line indicates the phase space background 
contribution. (b)The Dalitz plot. 
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Figure 5. (a). The m K -^ invariant mass spectrum 
from J/tp —* pK~K. (b). The m K -^ — m K - — 
mj^ distribution after the efficiency and phase space 
correction. 



ciency and phase space correction presents a 
more obvious peak at the K~ A mass thresh- 
old. 

The Partial Wave Analysis (PWA) is ap- 
plied to J/tp —> pK~A to study this near 
threshold enhancement, denoted as N x . Dif- 
ferent combinations of N* and A* states, 
listed in PDG, are included in m K -x and 
rripK- for the PWA fit. The results show 
that the mass and width of the enhancement 
are around 1500 - 1650 MeV and 70 - 110 
MeV, respectively, the J p of l/2~ is favored 
and the product branching ratio Br( J/ip — > 
pN x ) x (N x -> K~A) is around 2 x 10~ 4 . The 
big product branching ratio indicates a large 
coupling of N x to KA. 



4 Summary 

Based on 5.8 x 10 7 J ftp and 1.4 x 10 7 ip(2S) 
events accumulated at the BESII detector, 
no 9(1540) is observed in both J/tp and 
ip(2S) hadronic decays. From a sample of 
5.8 x 10 J/ip events, a narrow enhancement 
near 2m p in the invariant mass spectrum of 
-> jpp decays is observed, 
c.c. decays, an enhancc- 
u p + M\ mass threshold is 
observed in the combined pA and pA invari- 
ant mass spectrum. We also observed an en- 
hancement near the uik+Ma mass threshold 
in the same channel. The J p of this enham- 
cement being l/2~ is favored and it has a 
large coupling to KA. 



pp pairs from J /ip 
In J/tp -> pK-A - 
mcnt near the m n 
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